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In these lectures I highlight some key features of massive neutrinos in the context of 
cosmology. I first review the thermal history and the free-streaming kinematics of 
the uniform cosmic background neutrinos. I then describe how fluctuations in the 
phase space distributions of neutrinos and other particles arise and evolve after 
neutrino decoupling according to the linear perturbation theory of gravitational 
instability. The different clustering properties of massive neutrinos (aka hot dark 
matter) and cold dark matter are contrasted. The last part discusses the nonlinear 
stage of gravitational clustering and highlights the effects of massive neutrinos on 
the formation of cosmological structure. 



1 Neutrino Masses 

These lectures discuss how the universe serves as a learning ground for massive 
neutrinos. Before doing so, let us briefly review some experimental measure- 
ments of neutrino masses. 

Upper bounds on neutrino masses from kinematic measurements in labo- 
ratories continue to improve. El For the r-neutrino, m^^ < 18.2 MeV from the 
decay channel r — > Stt + i>r. For the /x-neutrino, m^^ < 170 keV from two- 
body pion decay. For the electron neutrino, the quantity is measured in 
tritium beta decay by fitting the shape of the energy spectrum near the end- 
point. Experiments thus far have yielded nonphysical negative values for rrif, , 
indicating unexplained systematic effects in the measurements. A conservative 
upper bound is put at m^^ « 15 eV. The spread in arrival times of neutrinos 
from supernova explosions provides an independent way to constrain the mass 
of the electron neutrino. Various limita^have been reported for SN 1987A; a 
conservative estimate is < 23 eV.&El 



2 Properties of Cosmic Background Neutrinos 

2.1 Temperature and Density 

For a brief 1 second after the big bang, neutrinos enjoy being part of the 
thermal bath composed of photons, electrons, protons, neutrons, and the as- 
sociated anti-particles (after the quark-hadron era) . The weak interactions at 
this early time are rapid enough to keep these particles in thermal equilib- 
rium at a single temperature T. After 1 second, when T drops below about 1 
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MeV, however, the neutrino interaction rate becomes slower than the Hubble 
expansion, and neutrinos become effectively coUisionless and freely-streaming 
particles whose trajectories are determined by the geodesic equations. This 
event is commonly referred to as "neutrino decoupling." As the universe ex- 
pands, the momenta and temperature of neutrinos are simply redshifted, and 
the neutrino temperature is given by the familiar formulas 

T,(a) = a- ir,,o, ^-."=(yt) T-,^o^lM7K, (1) 

where a is the cosmic scale factor, the subscripts denote the present-day 
values, aijd the cosmic background photon temperature is taken to be Tj_o — 
2.728 KM 

An important feature of the neutrino distribution after decoupling is that, 
although weak interactions are no longer rapid enough to keep neutrinos in 
thermal equilibrium with other particle species, neutrinos retain their equi- 
librium distribution as long as no other physical processes (e.g., gravitational 
clustering; see Sec. 3) are present to alter it. Therefore, to zeroth order in den- 
sity and metric perturbations, the phase space distribution /o of the cosmic 
background neutrinos is of the simple Fermi-Dirac form 

where e = a(p^ -I- m^)^/^ is the comoving energy, T^^ is the neutrino tempera- 
ture given by Eq. (1), is the number of spin degrees of freedom, and hp and 
ks are the Planck and the Boltzmann constants. 

The situation is further simplified if neutrino masses are ^ 1 MeV. Such 
neutrinos are highly relativistic at decoupling; their energy e, and hence the 
distribution function /o, are independent of to a good approximation. One 
can easily show that, as long as rrii, ^ 1 MeV, the number density of the 
cosmic background neutrinos is related to the neutrino temperature by 

where C(3) ~ 1.202 is the Riemann zeta function of order 3. This gives a 
present-day density of w 113 cm~'^ for every neutrino species independent of 
their masses. (For comparison, the present-day photon density is « 412 cm~^.) 
It also follows that the contribution of these neutrinos to the present-day mass 
density parameter, O^, is related to their masses by the simple relation 



where the index i runs over all light, stable neutrino species (e.g., Vf.^v^^ and 
Vr), and the Hubble constant is Hq = 100 ft. km s^^ Mpc^^. One then arrives 
at the important conclusion that in order for neutrinos not to close universe 
(i.e. Vly < 1), the sum of neutrino masses must not exceed 93 eV. This value 
is far below the current laboratory limits (see Sec. 1). Cowsik & McClelland □ 
were the first to use such cosmological arguments to place an upper bound on 
neutrino masses. (Unfortunately, these "hot dark matter" models in which the 
mass density is dominated by massive neutrinos have been found to produce 
excessive large voids surrounded by large coherent sheets and filaments that 
are not seen in the observable universe, cl Modifications to this model will be 
discussed below.) 

In the high mass regime, rrii, 3> 1 MeV, there exists another window where 
the neutrino contribution to the mass density parameter f2 of the universe is 
subcritical. The argument is that neutrinos with rrii, 3> 1 MeV become non- 
relativistic long before decoupling. Neutrino and anti-neutrino pairs cease to 
be created in abundance once the thermal temperature drops below m^^ and 
the neutrino density is suppressed by the Boltzmann factor e"™'^/'^^ . This 
large reduction factor in the relic abundance allows neutrinos to hajffi large 
masses without overdosing the universe. A more careful calculation eI shows 
that an < 1 universe implies a lower limit of ~ 2 GeV if these heavy neutrinos 
are Dirac, and ~ 6 GeV if they are Majorana. Since this mass range is well 
above the current upper mass bounds from laboratory measuremenls, it is of 
interest only when one considers more exotic theories for neutrinos. Q 



2.2 Kinematics and Free Streaming 

Let us now turn to the kinematics and the streaming properties of neutrinos. 
In general, neutrinos of mass rui, become non-relativistic after a redshift of 

This redshift has important implications for structure formation because it 
dictates the time at which massive neutrinos begin to make a transition from 
being radiation to matter. Note that this transition occurs fairly early, before 
recombination if rrijy ^ 1 eV. The average momentum of the cosmic background 
neutrinos at temperature Ti, is given by 



{p) = 3.15kBT„/c. 
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(6) 



In the non-relativistic regime (p = mi,v), the average neutrino speed can be 
written as 

Since Tj, cx a^^ cx (1 + z), massive neutrinos slow down as time goes on. It 
is important to keep in mind that neutrinos with a mass of several eV have 
slowed down to an average velocity below 100 km s~^ today. 

We also note that at the redshift of matter-radiation equality, Zcq ~ 
24000 f^/i^ (i.e. when the total energy density in radiation in the universe 
equals that in matter), light neutrinos with 1 < < 10 eV are zooming 
around with speeds close to c. Such large thermal speeds prevent massive neu- 
trinos from clustering gravitationally during this epoch, and this is why light 
neutrinos are referred to as hot dark matter (HDM). In contrast, perturbations 
in cold dark matter (CDM), which by definition has negligible thermal veloc- 
ities, can grow unimpeded after Zoq. I will quantify the different clustering 
behavior of CDM and HDM further in Sec. 3 and 4. 

Since neutrinos cannot cluster appreciably via gravitational instabilities on 
scales below the free streaming distance, this introduces a characteristic length 
scale into the problem. This scale is given by the free-streaming wavenumber 
(in comoving coordinates) 

4 = ^, (8) 

{v) 

which is analogous to the Jeans length for a self-gravitating system of den- 
sity p. For k < kfs (i.e. large wavelengths) , the density perturbation in the 
neutrinos is Jeans unstable and grows unimpeded in the matter-dominated 
era. For k > kfs, the density perturbation decays due to neutrino phase mix- 
ing. A phase-space interpretation of the free streaming property is that the 
phase Jnixing of coUisionless particles damps the growth of density perturba- 
tions. 19 When the neutrinos are relativistic, (w) « c, and the free-streaming 
distance is approximately the particle horizon, which scales as /cts(a) cx (in 
the radiation-dominated era). After the neutrinos become non-reJativistic the 
relations (v) cx a~^m~^, rrii, cx fli^h'^, and p cx a^^h'^ then implyEI 

/cfe(a) cx a}l'^VL^h^ . (9) 

As expected, the free-streaming distance (cx fcg^^) decreases with time as the 
neutrinos slow down. This also implies that neutrinos, can cluster gravitation- 
ally on increasingly small length scales at later times. la Such behavior has been 
seen in cosmological numerical simulations and will be discussed in Sec. 4. 
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3 Linear Perturbations in Neutrinos and Other Particles 

Thus far our discussion has focused on the properties of the smooth cosmic 
background neutrinos, and Eqs. (l)-(7) were derived under this assumption. 
The universe today, however, is clearly far from being homogeneous on scales 
of ^ 100 Mpc and below. Baryons and dark matter in galaxies, clusters, and 
superclusters show a wide spectrum of overdensities above the cosmic mean. 
The current theoretical framework for the origin and evolution of these cos- 
mic structures rests upon the assumption that certain primordial fluctuations 
(perhaps originated from quantum fluctuations of scalar fields during the in- 
flationary era) imprint a perturbation spectrum on all matter and radiation. 
These fluctuations subsequently grow via gravitational instabilities to give rise 
to the wide range of observed structures. How are the cosmic relic neutrinos 
affected by all this? 

To understand the growth of density perturbations in neutrinos as well as 
other forms of matter and radiation, one would need to learn the linear cos- 
mological perturbation theory of gravitational instability. A full description 
of this theory requires more time than is allocated for these lectures. I will 
only sketch the theory below with emphasis on the neutrino component. In- 
terested readers should refet-to the pioneering work of Lifshitz, later reviewed 
in Lifshitz & Khalatnikov. E3 More modern treatments of various aanects of 
this theory can be found in the textbooks by Weinberg and Peebles, Eil in th£ 
reviews by Kodama & Sasaki and Mukhanov, Feldman & Brandenberger, — 
and in the Summer School lectures by Efstathiou, Bertschinger, and Bond.E^ 
A complete description of this theory for all relevant particles is given by Ma & 
Bertschinger. Eil Here, I will only discuss the essence of the theory and highlight 
the physical meaning of the key results. 

3.1 Neutrino Phase Space 

Let us start with the neutrinos. The full phase space distribution function of 
neutrinos can be written as 

/(f,p,i) = /o(p) + /i(^,P,i), (10) 

where /i denotes perturbations to the Fermi-Dirac distribution /g given by 
Eq. (^. Unlike the unperturbed term /q that depends only on p, fi can 
have complicated dependence on time as well as positions x and the conjugate 
momenta p. The equations for neutrino temperature, number density etc. 
discussed in Sec. 2 were obtained assuming f — fo- A non- vanishing /i would 
lead to perturbations in these quantities. For example, the perturbed neutrino 
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energy density is related to /i by 

Sp{x,t) = / d^pe fi{x,p,t) 



(11) 



Other quantities such as perturbations in the pressure and shear can also be 
related to fi. 

It is in general difficult to compute and sample /i directly because at 
a given time, it depends on six variables. A Monte Carlo technique, or a 
"general-relativistic A^-body" technique, has been developed to evolve fi from 
redshift z ^ 10^ shortly after neutrino decoupling until z ~ 10 when nonlinear 
effects become non-negligible. cB In this calculation, an ensemble of neutrino 
simulation particles is initially assigned velocities drawn from the Fermi-Dirac 
distribution, which is an excellent approximation at z ~ 10^. The trajectory of 
each neutrino simulation particle is then followed by integrating the geodesic 
equations in the perturbed background spacetime. The metric perturbation 
gives rise to a nonzero /i, and the particle positions and velocities at a later 
time t represent a realization of /i (x, p, t) . Results from this calculation have 
revealed that at z ~ 15, positive correlations have developed in the rms neu- 
trino velocities and the overdensity, which would be absent if the phase-space 
distribution were purely Fermi-Dirac (i.e. / — /o). The more spatially clus- 
tered neutrinos are found to move faster, possibly resulting from an increase 
in the kinetic energy during gravitational infalls. 

3.2 Evolution of Perturbed Density Fields 

The phase-space description above is applicable to all particle species and is 
used in the full theory. The full, general-relativistic version of the linear per- 
turbation theory is described by a set of coupled and linearized Einstein, Boltz- 
mann, and fluid equations. The variables include the metric perturbations to 
the homogeneous and isotropic Friedmann-Robertson- Walker spacetime, and 
the phase-space perturbations in all relevant particle species (e.g., photons, 
baryons, cold dark matter, massless and massive neutrinos). The Einstein 
equations describe how the time evolution of the metric perturbations is af- 
fected by the perturbations in the density, pressure, shear, and higher-order 
moments of matter and radiation. The Boltzmann and fluid equations, on the 
other hand, describe the time evolution of the radiation and matter distribu- 
tion in the perturbed spacetime. Together, this theory describes the growth of 
metric and density perturbations throughout the early history of the universe, 
and it serves as the foundation for all calculations of the linear power spec- 
tra for matter and temperature variations imprinted on the cosmic microwave 
background. 
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Figure 1: Time evolution of the perturbed energy density field, 5 = p/ p — 1, for five mat- 
ter and radiation components in a flat C+HDM cosmological model. (See text for model 
parameters.) The results are from integration of the coupled Einstein and Boltzmann equa- 
tions. Since the equations are linearized, each fc-modc evolves independently. Two modes 
are shown here for illustration. In each panel, the five curves represent 5 for the cold dark 
matter (solid), baryons (dash-dotted), photons (long-dashed), massless neutrinos (dotted), 
and massive neutrinos (short-dashed), respectively. 
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Figure illustrates a small subset of results that can be obtained from 
numerical integration of these linearized equations. It shows the time evolution 
of the perturbed energy density field, 6 — 6p/p = p/p— 1, for the five relevant 
particle species in a cold+hot dark matter (C+HDM) model. This model 
assumes an Einstein-de Sitter universe containing a mixture of CDM, HDM (i.e. 
massive neutrinos), baryons, photons, and massless neutrinos. The density 
parameters in the first three components in this model are ilc = 0.75, il^ = 
0.2, and fib — 0.05, and the Hubble parameter is taken to be Hq = 50 km 
s~^ Mpc~^, or ft, = 0.5. From Eq. (4), these parameters correspond to a 
neutrino mass of 4.7 eV. (For definiteness, this calculation has assumed that 
only one type of neutrinos, presumably Vt, has a non-negligible mass.) Two 
wavenumbers, k = 0.01 Mpc~^ (top) and 1 Mpc~^ (bottom), are shown in 
Figure |l| to demonstrate the intricate dependence of 5 on length and time 
scales. The overall normalization of 5 is set arbitrarily. 

We observe several salient features in Figure |l|. First, the amplitudes of 5 
for all particles grow monotonically until a critical time, after which different 
particle species exhibit very different behavior. This critical time is the "hori- 
zon crossing" time, and it occurs when the horizon has grown large enough to 
encompass the wavelength of a given mode of perturbation. A mode of per- 
turbation is therefore not in causal contact until horizon crossing. Naturally, 
this occurs earlier for smaller wavelengths (i.e. larger k). In Figure |l|, one 
can indeed see that the A: = 1 Mpc~^ mode enters the horizon at a ~ 10~^ 
while the k = 0.01 Mpc~^ mode enters the horizon a little after a = 10"'^. 
A point to keep in mind is that the behavior of 5 before horizon crossing is 
strongly dependent on the choice of gauge. The results shown in Figure |l] are 
computed in the so-called synchronous gauge, which is a popular choipe due 
to historical precedent. t£l See Mukhanov, Feldman & Brandenberger t3 and 
Ma & Bertschingerliil for discussion of a more convenient gauge (the conformal 
Newtonian gauge). 

The second feature in Figure to note is that after horizon crossing, the 
photons (long-dashed) and baryons (dot-dashed) exhibit rapid, coupled oscil- 
lations in the k — \ Mpc^^ mode but only the photons oscillate in the k = 0.01 
Mpc~^ mode. This occurs because the former enters the horizon before recom- 
bination at flj-cc ~ 10"'^, and the photons and baryons are coupled by Thomson 
scattering and oscillate acoustically. The coupling is not perfect. The_friction 
of the photons dragging against the baryons leads to Silk damping, E3 which 
is prominent in the bottom panel of Figure ^ at a ~ lO""^'^. After recombina- 
tion, the baryons decouple from the photons and fall quickly into the potential 
wells formed earlier by the CDM. This results in the rapid growth of the dot- 
dashed curve in the bottom panel of Figures 111 The mode with k = 0.01 
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Mpc^^, on the other hand, enters the horizon when the universe is neutral. 
The baryons therefore grow hke the CDM and do not osciUate. The critical 
length scale demarcating these two regimes is the horizon size at recombination 
Oioc: fcioc ^ 0.03 Mpc"^ 

The third feature to note in Figure |^ is the rate of growth of the CDM 
component (solid curve) after horizon crossing. Close inspection shows that 
the CDM in the bottom panel grows more slowly at 10^''' ^ a <^ 10^^ than later 
on, whereas in the upper panel, the CDM simply grows with a power law after 
horizon crossing. This is because the shorter wavelength mode (fc = 1 Mpc~^) 
enters the horizon when the energy density of the universe is dominated by ra- 
diation, and fluctuations in matter (e.g. CDM) cannot grow appreciably during 
this era. The critical scale separating continual and suppressed growth is the 
horizon size at the time of radiation-matter equality Ocq ^ 4 x 10^^ (ilh^)^^ : 
fceq ~ 0.1 Mpc"'^ for the parameters of this model. 

The fourth feature to note in Figure |l] is the behavior of the massive 
neutrinos (short-dashed). As discussed in Sec. 2, neutrinos of masses within 
the cosmologically interesting range (~ 1 to 10 eV) are highly non-relativistic 
today but were relativistic at earlier times. This property is in fact evident 
in the top panel of Figure |l|: Careful inspection shows that at a « lO"'', the 
density field S in massive neutrinos is indeed making a gradual transition from 
the upper line for the radiation fields to the lower line for the matter fields. 
(More precisely, the primordial perturbations are assumed to be "isentropic" 
here, which leads to a perturbation amplitude that is a factor of 4/3 higher for 
radiation than matter.) The subsequent evolution of S in massive neutrinos for 
this mode (fc ~ 0.01 Mpc^^) is very similar to that of CDM. This is because it 
enters the horizon when the thermal velocities of the neutrinos have decreased 
substantially; the free-streaming effect is therefore unimportant. For the k = 1 
Mpc~^ mode, on the other hand, the free streaming effect is evident and the 
growth of S in massive neutrinos is suppressed until the characteristic free- 
streaming scale A;fs(a) given by Eq. (9) grows to ~ k. Afterwards, the short- 
dashed curve for massive neutrinos is seen to grow again and catch up to the 
CDM. Since fcts oc a^/^, the larger k modes suffer more free-streaming damping 
and 5 for massive neutrinos can not grow until later times. The damping in 
the massive neutrino component also affects the growth of the CDM, slowing 
it down more for models with larger compared to the pure CDM model. 

3.3 Growth Rate and Power Spectrum 

I have used Figure ^ computed for two particular fc-modes in a particular 
cosmological model to illustrate the physical meaning of many key features 
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in the evolution of the density field for matter and radiation throughout the 
cosmic history. I will now discuss general descriptions that can be conveniently 
used to characterize the fluctuation amplitudes over a range of length scales 
for a variety of cosmological models. For example, it is extremely useful to 
know the dependence of S and its time derivative on A; at a given time for 
a wide range of models. The most basic quantity to use is the linear power 
spectrum, P{k,t), and the growth rate of the density field, / = (ilog(5/(iloga. 
{Note: I am following the convention of using f to denote the growth rate; it 
should not be confused with the phase-space distribution function of Sec. 2 and 
3.1.) The power spectrum quantifies the two-point statistics of 5, and for a 
Gaussian field, P{k) represents its rms fiuctuations and completely specifies 
its statistical properties. The power spectrum is therefore of fundamental 
importance in cosmology. 

Comparing the growth rate and power spectrum for models with and with- 
out massive neutrinos is an effective way to illustrate the effects of hot dark 
matter. In the standard CDM model with = 1 and h = 0.5 (neutrino mass 
is assumed to be zero), the CDM density field grows as the expansion factor 
a on all scales; therefore / = 1. As discussed in Sec. 2.2, massive neutrinos 
introduce an additional length scale, the free-streaming distance, below which 
fiuctuations are washed out and the growth rate is retarded. The growth rate 
is therefore more intricate in models with massive neutrinos and is generally 
a function of the wavenumber fc, neutrino density parameter fi^, and time. 
Figure || illustrates such dependence in four different C+HDM models that 
assume a mixture of CDM and HDM. It shows that the growth is suppressed 
at large fc, and models with a larger fraction of energy density in HDM suffer 
more. It also shows that the suppression becomes less severe at later times. 

One can gain some understanding of the behavior shown in Figure ^ by 
exploring two asymptotic regions that can be solved analytically for C-I-HDM 
models with ^Ic + — (1) Since HDM behaves like CDM above the free- 
streaming distance, / — + 1 as fc — + 0; (2) In the opposite limit of large fc, the 
HDM density field 6h is severely dampened compared to the CDM density field 
5c because of the neutrino effects. For 5h <C 5c, the time evolution of the CDM 
density field is governed by the linearized fluid equation 

5c + -5c = l.bH^a^nc5c , (12) 
a 

where the dots denote differentiation with respect to the conformal time r. 
Since Ha = 2/t in the matter-dominated era, the growing solution in this 
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Figure 2: Growth rate of the CDM density field, / = dlog <5/(i log a, in four flat C+HDM 
models at a = 1 (solid) and 0.1 (dashed). The four models assume different neutrino masses: 
m^y = 1.2, 2.3, 4.6, and 6.9 eV (from top down), corresponding to f2,y = 0.05, 0.1, 0.2, and 
0.3. At small fc, the CDM density field in these models grows with the same rate {5 cc a) as 
in the standard CDM model. At large k, the growth rate is suppressed because a fraction of 
the energy density in C+HDM models is in the hot neutrinos that exhibit less gravitational 
clustering. The suppression at large k becomes less severe at later times because the velocities 
of the hot neutrinos decrease with time. 
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regime is easily shown to be 







^ fik ^ oo) = + 24f7, -^ = yi-^n,-l. (13) 
It is interesting to note that 

/oo « f^:?-' (14) 

is an excellent approximation to the equation above, especially for the cosmo- 
logically viable range of il^, <, 0.3. Using these analytic solutions in the asymp- 
totic regimes and the scaling dependence of the free streaming wavenumber kfg 
in Eq. (^, one can construct a simple approximation for / for a wide range 
of model parameters. It is found that the growth rate / is well approximated 
byB 

rflog^ _ l + f^°-6 0.00681x1-363 _ k 

■' " dhi^ ~ 1 + 0.00681x1-363 ' = IVI ' ^ ' 

where V^, is a shape parameter derived from Eq. (^, 

= a^l'^VL^h^ , (16) 

fic + iliy = 1, and k is in units of Mpc^i. Note that Eq. ( p^ depends only 
on the variable x that characterizes the neutrino free-streaming scale, and 
(or fij,) via /oo- The fractional error of the fit relative to the numerically 
computed values is smaller than 0.5% for a wide range of parameters. The 
seemingly complicated multi-parameter dependence of Figure ^ is succinctly 
incorporated in Eq. (p^). 

For the linear power spectrum, the slower time growth of 5 at fc > k^s 
for the C-t-HDM models shown in Figure ^ indicates a suppressed clustering 
amplitude on these scales. This effect is illustrated in Figure pi which shows 
the density-averaged power spectrum (i.e. P = {^c^fPc + ^v^/Pv]^) for three 
flat C-t-HDM models. In general, Pu ^ Pc on small length scales due to the 
neutrino thermal velocities, and the models with higher clearly have less 
power at large k in accordance with Figure |[ A good approximation for the 
relative P(fc) in a flat^C-hHDM model (with Q.^ < 0.3) and a pure CDM model 
{fli, = 0) is given hyU 

P{k, a, n^) (1 + bi J'^l'^ + 62 \ " 



p{k,a,n, = o) y l + b^xl" J 

k 

x = — , a;o = a;(a=l), (17) 

-L ij 
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Figure 3: Power spectrum at z = 1.5 for the perturbed matter density field for three flat 
C+HDM models with = 0.1, 0.2, and 0.3. The solid curves show the linear P(fc) com- 
puted from the linear perturbation theory; accurate analytical fitting formulas are given 
by Eq. (p^. The triangles show the nonlinear P(fc) computed from Af-body simulations; 
accurate analytical approximations are given by Eq. (|lq). 
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where the best-fit parameters are bi = 0.004321, 62 = 2.217 x lO^^^ 5^ ^ 11.63, 
and 64 = 3.317 for k in units of Mpc^^. Analytical approximations for the 
separate cold and hot spectra Pc and can be found in the same reference. 
More complicated aiiproximations for a wider range of parameter space have 
also been proposed. Il3 

4 Nonlinear Gravitational Clustering of Neutrinos 

As far as cosmological structure formation is concerned, the main difference 
between massive and massless neutrinos is that the former can participate in 
the processes of gravitational clustering and hence serves as a component of 
the dark matter in the universe. Massless neutrinos, on the other hand, affect 
cosmology only through their contribution to the radiation energy density. Any 
primordial perturbations in massless neutrinos are damped out after horizon 
crossing as a result of phase mixing (see dotted curves in Figure |l|) , and the 
only remnant of this component is the elusive T^^q — 1.947 K background 
described in Sec. 2. 

^.1 Spatial Distribution of Neutrinos 

For massive neutrinos in models with a mixture of CDM and HDM, it is 
interesting to ask: do the neutrinos fall in the CDM potential wells and form 
a part of dark matter halos? One may naively think not because neutrinos are 
too hot. Our discussion thus far, however, indicates otherwise. We have seen 
in Eq. (^) that neutrinos slow down with the expansion of the universe. Those 
with a mass of several eV are travelling with a speed much below the typical 
velocity dispersions ^ 200 km s^^ of stars in galaxies. They can potentially 
be bound to galactic halos. In general, the extent to which massive neutrinos 
can cluster gravitationally depends on their mass and speed. 

To illustrate this point further, let us examine results from numerical sim- 
ulations. Figure ^ shows the projected spatial distribution of cold particles 
(middle), hot particles (bottom), and the sum of the two (top) in a simulated 
dark matter halo in a flat fli, = 0.2 C+HDM model (m,, = 4.7 eV). Three red- 
shifts are shown (from left to right): 2 = 2,1 and when the universe is 2.5, 4.6, 
and 13 billion years old, respectively. Each panel is 3.5 x 3.5 h^^ Mpc in phys- 
ical coordinates. The parent simulation is a large iV-lipdy run with 23 million 
simulation particles in a (100 Mpc)'^ comoving box. Ej The dense halo shown 
at 2: = is clearly formed from mergers of two smaller halos and their satel- 
lites at higher redshifts, demonstrating the "bottom- up" hierarchical pattern 
of structure formation which is preserved in C-I-HDM models with flc > iljy. 
Massive neutrinos are visibly clustered in the bottom panels, but their spatial 



14 



Figure 4: Spatial distribution of cold dark matter (middle) vs. hot dark matter (bottom) 
particles in a simulated halo formed in a large A'^— body run for the flat Q,y = 0.2 C+HDM 
model. The top panels shows the sum of cold and hot particles. Three rcdshifts arc shown: 
z = 2, I, and (from left to right), corresponding to cosmic times of 2.5, 4.6, and 13 billion 
years. All boxes have the same physical scale, 3.5 x 3.5 Mpc. The final halo at z = is 
clearly a merger product of two dominant subhalos and several smaller satellites at higher 
redshifts. The spatial distribution of the massive neutrinos (4.7 eV) is visibly smoother than 
that of the CDM. 
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distribution is smoother than that of the CDM. In the smaller CDM clumps at 
z ^ 1, there are no discernible HDM halos at the same locations. The average 
thermal speed of an ensemble of 4.7 eV Fermi-Dirac neutrinos at z ^ 1 is ^ 70 
km/s (cf. Eq. |^). As discussed in Sec. 3.1, when linear perturbations are taken 
into account, neutrinos in overdense regions have even higher velocities. It is 
therefore not surprising that the shallower potential wells of these small halos 
cannot trap a substantial number of HDM particles. 

The clustering of the 4.7 eV neutrinos shown in Figure ^ should be com- 
pared with the constraint derived by Tremaine & Gunn,LJ which states that 
if cosmic neutrinos were to make up the bulk of galactic and cluster halos, 
they must be more massive than ~ 10 eV so that the Pauli exclusion principle 
is not violated. The C+HDM models considered here assume ilc > so 
as to preserve the successful hierarchical formation of structure in pure CDM 
models. Dark matter halos such as in Figure ^ have a substantial fraction of 
CDM, which in turn enhances the gravitational infalls of massive neutrinos. 
The clustering of neutrinos in C+HDM models is therefore more complicated. 



4-2 Nonlinear Power Spectrum 

The process of nonlinear gravitational clustering can be quantified statistically. 
Here I will only discuss the lowest-order statistical description given by the 
nonlinear power spectrum P(fc) of the matter density field. 

We have already discussed the linear power spectrum in Sec. 3.3. The 
triangles in Figure ^ show the nonlinear P(fc) computed from the particle 
positions in numerical simulations of three C-I-HDM models. The hierarchical 
nature of gravitational collapse in these models is illustrated by the fact that 
the high-fc modes have become strongly nonlinear whereas the low-fc modes 
are still following the linear power spectrum. The fact that the lowest several 
k modes are still linear at z = ensures that our choice of the simulation box 
size (100 Mpc) is large enough to include all waves that have become nonlinear 
at present. 

The calculation of the fully evolved P{k) for a given model is a laborious 
task involving the execution of high-resolution A^-body simulations. Fortu- 
nately, some recent progress has been made in constructing analytical fitting 
formulas for a wide range of interesting models. The strategy is to examine the 
mapping between the linear and nonlinear P{k) for a small, selected set of mod- 
els with A^-body data and then to extract systematic behavior for a wider range 
of parameters. The work carried out,th.us far has investigated scale- free models 
with a power-law power spectruni,K4 ^pure CDM and CDM wiith a cosmolog- 
ical constant A (LCDM) models, BNel and C-^HDM models.EJ The proposed 
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formulas typically have a functional form that is motivated by analytical solu- 
tions in asymptotic regimes, but in order to obtain accurate approximations, 
the coefficients are calculated from fits to the nonlinear P{k) computed from 
the numerical simulations. These approximations have provided physical in- 
sight into the process of nonlinear collapse and much practical convenience in 
incorporating the prominent nonlinear effects illustrated in Figure ^ 

The formula applicable for the widest range of cosmoiogicai models thus 



far is given below. It mapj 
and nonlinear regimes byr^ 



he density variance A(fc) = Ank P(k) in the linear 



Ani{k) Ai{ko) 



G x = 1 + ln l + 0.5a; — ^\ . 18 



Note that Ai and Ani are evaluated at different wavenumbers, where /cq — 
A: (1 -I- Aui)^^/"^ corresponds to the pre-coUapsed length scale of k. The pa- 
rameter (jg is the rms linear mass fluctuation on 8h~^ Mpc scale at the red- 
shift of interest, and (3 = 0.7 + IQflf,. The functions go = g{il.^,fl\) and 
g = g(f2in(a), ^tSa)) are, respectively, the relative growth factor for the linear 
density field at present day and at a for a model with a present-day matter 
deusity fim and a cosmoiogicai constant IIa- A good approximation is given 



byB 

g = ^l]„,(a)[r!,„(a)4/7 - ^^{a) + (1 + an(a)/2) (1 + r!A(a)/70)]-i , (19) 

and f]m(a) = f^m a" V[l + ^m(a"^ - 1) + f^A(a^ - 1)] and ^tSa) =nAaV[l + 
f^m(a~^ — 1) + ^K{a^ — 1)]- The time dependence is in factors (Xg and g. 
For CDM and LCDM models, a good fit is given by ci = 1.08 x 10"'* and 
C2 = 2.10 X 10"^ For C+HDM, a good fit is given by Ci ^ 3.16 x 10^^ ^nd 
C2 = 3.49 X 10^4 for 9.^ = 0.1, and ci 6.96 x lO'^ and ca = 4.39 x lO"'^ for 
= 0.2. 



^.3 High-Redshift Constraints 

The C+HDM models discussed thus far are a class of models bridging the 
much-studied albeit troubled pure CDM and pure HDM models. They are 
parameterized by the neutrino density parameter < rj^^ < 1, or equivalently, 
by the neutrino mass < < 93/i^ eV (see Eq. (4)). The original motivation 
for examining these mixed models is to study whether the free-streaming effect 
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introduced by the massive neutrinos could suppress the growth of density per- 
turbations below the free-streaming scale, and thereby alleviate the problem 
of excess small-scale clustering in the standard CDM model. As illustrated in 
Figure ||, massive neutrinos do indeed reduce the amplitude of clustering at 
large fc, and a larger fi^ leads to smaller high— fc power. 

Any viable cosmological model that provides a good statistical match to 
the local universe must also reproduce the appropriate evolutionary history 
out to high redshifts. Although gravitational clustering on galactic scales is 
indeed reducfidJn C-f HDM models, providing a better match to low-redshift 
observations pEj this suppression in the clustering power at high redshifts has 
been found to pose serious problems for some models. Studies based on semi-, 
analytic theories and dissipationless numerical simulations have shown tSEZi 
that flat C-I-HDM models with ^l^ > 0.2 do not produce enough early structure 
to explain the statistics of damped Lya systems at redshift z > 2. More 
recent workcl has included the effects of gas ionization and dissipation in the 
theoretical calculations and has compared the results to new data for damped 
Lya systems at even higher redshift z ~ 4. It is found that ionization of 
hydrogen in the outskirts of halos and gaseous dissipation near the halo centers 
tend to exacerbate the problem of late galaxy formation. The amount of dense 
gas associated with the damped systems falls well below that observed, even 
for the flat fijy = 0.2 C+HDM model. This has placed an upper bound of ^ 5 
eV on the sum of z/g, t'^, and Vr masses, which is much more stringent than 
the upper limits given by current particle experiments (see Sec. 1). 

Let us briefly discuss the limitations and uncertainties in these calculations. 
Still debated is the nature of damped Lya absorption - whether it is due 
to intervening large, rapidly-rotating disk galaxies c2l with circular velocities 
200 km s~^, or iofalling and randomly moving protogalactic gas clumps in 
dark matter halos t£l with virial velocities of ^ 100 km s~^. When kinematic 
considerations are included, the latter model may have pmblems balancing 
the high energy dissipation rate caused by cloud collisions. Eil This uncertainty 
aside, in either scenario, a host dark halo of velocity at least ^ 100 km s^^ is 
needed to reproduce the large velocity widths and asymmetries of the observed 
low-ionization lines associated with Lya systems. Uncertainties associated 
with the finite resolution of simulations have also been studied in some detail. o 
It is found that even when the contribution from the numerically unresolved 
halos with velocities v <J 100 km s^^ is included, the absorption incidence is 
increased by at most a factor of 2. This is insufficient to erase the discrepancies 
reported for flat C-t-HDM models with fl^, > 0.2. The upper bound on neutrino 
masses from current cosmological studies is therefore ~ 5 eV. 
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